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ABSTRACT 

We present a chemical abundance distribution study in 14 a, odd-Z, even-Z, light, and Fe- 
peak elements of approximately 3200 intermediate metallicity giant stars from the APOGEE 
survey. The main aim of our analysis is to explore the Galactic disk-halo transition region 
between -1.20 <[Ee/H] < -0.55 as a means to study chemical difference (and similarities) 
between these components. In this paper, we show that there is an a-poor and a-rich sequence 
within both the metal-poor and intermediate metallicity regions. Using the Galactic rest-frame 
radial velocity and spatial positions, we further separate our sample into the canonical Galactic 
components. We then studied the abundances ratios, of Mg, Ti, Si, Ca, O, S, Al, Ch-N, Na, 
Ni, Mn, V, and K for each of the components and found the following: (1) the a-poor halo 
subgroup is chemically distinct in the a-elements (particularly O, Mg, and S), Al, Ch-N, and 
Ni from the a-rich halo, consistent with the literature confirming the existence of an a-poor 
accreted halo population; (2) the canonical thick disk and halo are not chemically distinct in 
all elements indicating a smooth transition between the thick disk and halo; (3) a subsample of 
the a-poor stars at metallicities as low as [Ee/H] ^ -0.85 dex are chemically and dynamically 
consistent with the thin disk indicating that the thin disk may extend to lower metallicities 
than previously thought, and (4) that the location of the most metal-poor thin disk stars are 
consistent with a negative radial metallicity gradient. Einally, we used our analysis to suggest 
a new set of chemical abundance planes ([a/Ee], [Cn-N/Ee], [Al/Ee], and [Mg/Mn]) that may 
be able to chemically label the Galactic components in a clean and efficient way independent 
of kinematics. 


1 INTRODUCTION 

The stars in the Milky Way change their properties over time in 
almost every aspect. On one hand, their atmospheric parameters 
change as they age, producing significant changes in their colours 
and luminosities, with complete generations of stars enriching the 
interstellar medium with chemical elements as they die and new 
generations forming from this enriched material. On the other hand, 
their positions and kinematics change as they move with orbits gov¬ 
erned mostly by gravitation. The disentanglement of their spatial 
and dynamical properties, complemented with their ages and chem¬ 
ical compositions (when possible), are the tools that have been used 
for decades to classify the stars into different stellar populations 
building up a picture of the structure of the Milky Way we have 
today (e.g. Gilmore, Wyse & Kuijken 1989; Freeman & Bland- 
Hawthom 2002; Helmi 2008; Rix & Bovy 2013, and references 
therein). 

This classical picture of the Milky Way reveals a rather com¬ 
plex galaxy, with four primary “internal” components - the bulge, 
the thin disk, the thick disk and the halo, and one “external” com¬ 
ponent - the accreted material. These components form the basis 
of any attempt to model, to simulate and to understand the forma¬ 
tion and evolution of our Galaxy (see e.g. the reviews of Free¬ 
man & Bland-Hawthorn 2002; Ivezic, Beers & Juric 2012, and 
references therein). These components are usually defined by the 
phase spaces in which we observe them (spatially, kinematically. 


and chemically). However there is an ongoing debate as to whether 
these components are truly distinct objects or a part of a continuous 
sequence of an evolving galaxy (e.g. Bovy, Rix & Hogg 2012; Rix 
&Bovy 2013). 

Among the phase spaces to ‘separate’ the Galactic compo¬ 
nents, the method of using chemical abundance ratios as a means to 
probe the Galactic components has been performed and advanced 
significantly in the past decade by using relatively large samples 
of high resolution spectra (e.g. Edvardsson et al. 1993; Fuhrmann 
1998; Venn et al. 2004; Nissen & Schuster 2010; Sheffield et al. 
2012; Ramirez, Melendez & Chaname 2012; Feltzing & Chiba 
2013; Bensby, Feltzing & Oey 2014; Mikolaitis et al. 2014). In 
particular, the ratio of the mean of the a-elements^ relative to 
iron has been shown to be very powerful to disentangle the var¬ 
ious Milky Way components. Active research has been dedicated 
in order to explain the two sequences found in the abundances of 
the Q-elements respect to the abundances of iron. On one hand, at 
low metallicities ([Fe/H] ~ —1.2 dex), the 1ow-q sequence is at¬ 
tributed to stars that were chemically enriched in an environment 
separate from that of the bulk of the halo and later accreted on to 
the Galactic halo. At higher metallicities ([Fe/H] ~ —0.5 dex) the 
a-poor stars are interpreted as the thin disk. On the other hand, at 

^ e.g. Ca, Mg, Ti, O, S, and Si which are thought to be produced via Type 
II supemovae (SNII) 
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low metallicities the high-a sequence is attributed to the halo and 
at high metallicities the high-a sequence is attributed to the thick 
disk although there is likely significant overlap between the two. 

1.1 The canonical picture of the Galactic components 

Classically, the Galactic disk is often split into a thin disk and thick 
disk component (e.g. Gilmore & Reid 1983; Freeman & Bland- 
Hawthom 2002; Rix & Bovy 2013). The thin disk is thought to 
contain mostly young stars and is likely the final stage of dissipa¬ 
tive collapse. It is defined spatially by an exponential power law 
with a small vertical scale height and large radial scale length (e.g. 
Bovy et al. 2012; Haywood et al. 2013). Kinematically the thin disk 
stars follow near circular, co-rotational orbits with a low velocity 
dispersion (e.g. Edvardsson et al. 1993; Reddy et al. 2003; Kor- 
dopatis et al. 2013; Rix & Bovy 2013). Chemically, the thin disk 
its thought to extend over a metallicity range of -l-O.l < [Fe/H] < 
-0.70 dex and is near solar values in the a-elements (e.g. Bensby, 
Feltzing & Oey 2014). 

The thick disk was initially found to be separate from the thin 
disk on the basis of the spatial distribution of its stars (e.g. Yoshii 
1982; Gilmore & Reid 1983). It is thought to form potentially 
through a variety of processes including: satellite heating, accre¬ 
tion, or merging induced star formation, secular disk heating, radial 
migration, and others (for a review of these mechanism we refer the 
reader to Rix & Bovy 2013, and references therein). However the 
importance of each of these mechanisms in the formation and as¬ 
sembly of the thick disk is still not well understood. It is defined 
spatially by an exponential power law with a larger vertical scale 
height and smaller radial scale length compared to the thin disk 
(e.g. Bovy et al. 2012; Haywood et al. 2013). Kinematically, thick 
disk stars co-rotate with the disk, albeit with a smaller rotational 
velocity and overall have hotter orbits than their thin disk counter¬ 
parts (e.g. Edvardsson et al. 1993; Reddy et al. 2003; Bensby et al. 
2005; Haywood et al. 2013; Kordopatis et al. 2013). Chemically, it 
is thought to be distinct from the thin disk component in the [a/Ee]- 
metallicity plane with a high [a/Ee] signature (e.g. Bensby, Eeltz- 
ing & Oey 2014; Nidever et al. 2014; Recio-Blanco et al. 2014). 
The thick disk has a metallically that extends significantly lower 
than the thin disk, and can extend down to well below [Ee/H] < 
-1.0 dex (e.g. Beers et al. 2002). The thick disk is also thought to 
be older than the thin disk (e.g. Haywood et al. 2013; Masseron & 
Gilmore 2015). 

The Galactic halo is thought to naturally form very early on 
from a mixture of dissipative collapse and accretion (e.g. Eggen, 
Lynden-Bell & Sandage 1962; Searle & Zinn 1978; Ibata, Gilmore 
& Irwin 1994; Belokurov et al. 2006). As such most of its stars are 
very old (Jofre & Weiss 2011; Hawkins et al. 2014). It is spatially 
defined by a power-law with an index of approximately - 3.5 (e.g. 
Helmi 2008). Kinematically, the (inner) Galactic halo is thought 
to be pressure supported, have a small net prograde rotation (e.g. 
Carollo et al. 2010), and contain the highest velocity stars (e.g. 
Hawkins et al. 2015). Chemically, the Galactic halo is predomi¬ 
nantly metal-poor (e.g. Schlesinger et al. 2012) and enriched in the 
Q-elements (e.g. Ishigaki, Chiba & Aoki 2012). The (inner) Galac¬ 
tic halo has been shown to split into two chemically distinct com¬ 
ponents. Several recent studies have reported the presence of a- 
poor stars in halo samples (Nissen & Schuster 2010; Schuster et al. 
2012; Ramirez, Melendez & Chaname 2012; Sheffield et al. 2012; 
Bensby, Feltzing & Oey 2014; Jackson-Jones et al. 2014; Hawkins 
et al. 2014). These studies have shown that the a-poor sequence is 
distinct in kinematics, ages, and other chemical elements such as 


C, Na, and Ni compared to the a-rich sequence. It is thought that 
the a-poor sequence is assembled through the accretion of satellite 
galaxies. 

1.2 Decomposition of the canonical Galactic components 

Decomposing these components chemo-kinematically, particularly 
in the region of metallicity space where they overlap (-1.2 < [Fe/H] 
< -0.60 dex) is a significant challenge yet critical to understand 
the formation and assembly of the Galaxy, as this is the region of 
metallicity space where the thin disk, thick disk, and halo are all as¬ 
sembling and co-evolving. The problem arises with the fact that the 
original dynamical and spatial distributions of stars are perturbed 
over time, potentially erasing the memory of the original sites. 
Kinematical heating and spatial disruption can be produced via sev¬ 
eral process, such as bar resonances and radial mixing, dumpiness 
in the gas caused, and minor mergers, to name a few (e.g. Minchev 
et al. 2012; Rix & Bovy 2013; Haywood et al. 2013; Sanders & 
Binney 2013, 2015). Furthermore, dynamical and spatial distribu¬ 
tions depend heavily on the distance of the stars, which is subject of 
large uncertainties for the majority of them (e.g. Binney 2013, and 
references therein). Thus, matching models and simulations to ob¬ 
servational data continues to be one of the fundamental challenges 
in Galactic astronomy today. 

Despite this, the current mode through which Galactic com¬ 
ponents are decomposed is through kinematics via the Toomre di¬ 
agram with aid of the [a/Fe]-[Fe/H] plane. The left panel of Fig¬ 
ure 1 shows an illustration of the Toomre diagram which plots the 
quadrature of the vertical (denoted by W) and radial (U) veloci¬ 
ties as a function of the velocity along Galactic rotation (V). In 
this case, these velocities are relative to the local standard of rest 
(LSR). This diagram or a similar probabilistic kinematic Galac¬ 
tic component decomposition is widely used in the literature (e.g. 
Bensby, Feltzing & Lundstrom 2003; Venn et al. 2004; Nissen & 
Schuster 2010; Schuster et al. 2012; Ishigaki, Chiba & Aoki 2012; 
Bensby, Feltzing & Oey 2014; Hawkins et al. 2015). The power 
of the Toomre diagram lies in the fact that the canonical thin disk, 
thick disk, and halo have increasingly hotter kinematics (i.e. each 
component lies in an increasingly larger constant velocity circle on 
the Toomre diagram) and thus are further elevated in the diagram. 
The fundamental drawback to this diagram is the need for accurate 
proper motion and distances to fully resolve the 3D velocity vector. 

The [a/Fe]-meallicity plane (a cartoon version of this is 
shown in the right panel of Figure 1) is another approach to 
decompose the Galactic components. This method works because 
the only “easily measurable” property of long-lived stars that 
remains relatively unchanged over cosmic times is their chemical 
patterns that are imprinted in their spectra. Current large stellar 
spectroscopic surveys of the Milky Way, especially those of high 
resolution like APOGEE (Eisenstein et al. 2011), GALAH (De 
Silva et al. 2015) or Gaia-ESO (Gilmore et al. 2012), will enable 
methods such as chemical tagging stars providing us with a new 
picture for our Galaxy (e.g. Eeltzing & Chiba 2013). Under this 
scheme, a given stellar population is marked by the chemical 
distribution of the medium from which the stars formed. This 
distribution is principally defined by the initial mass function of 
the previous generation of stars that enriched the gas, the star 
formation rate and the different yields from the nuclear reactions. 

In this work, we study the stars at the Galactic disk(s)-halo 
interface by taking advantage of the APOGEE public spectro¬ 
scopic survey (Eisenstein et al. 2011), which contains not only 
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Figure 1. A) A cartoon illustration of the decomposition of the Galactic components along the Toomre diagram which plots the quadrature of the radial (U) and 
vertical (W) velocities as a function of the velocity along rotation (V). All velocities in the Toomre diagram are relative to the LSR. (B) A cartoon illustration 
of the decomposition of the Galactic components along the [a/Fe]-metallicity plane. 


thousands of stars at this interface, but its very high-resolution 
spectra allow us to go far beyond the study of only Q-elements. 
Specifically we are interested to know if the Galactic components 
must be defined in chemo-kinematic-spatial spaces or if they can 
be defined purely by chemistry. In this paper, we present a detailed 
chemical abundance evolution analysis of 14 different elements 
available in the APOGEE survey (including a, odd-even, light, 
and Fe-peak elements) with the aim to discuss how the thin disk, 
the thick disk, the halo and the accreted halo stars are linked. In 
particular, we are interested in the a-elements, carbon, nickel (and 
other Fe-peak elements), and other odd-even elements as they can 
be used to distinguish the accreted halo from the canonical halo, 
and canonical thin disk from the canonical thick disk (e.g. Bensby, 
Feltzing & Oey 2014). In section 2 we describe APOGEE data 
selection criteria and the subsample definitions. In section 3, we 
compare and contrast, element-by-element, the various Galactic 
components in order to better understand how they are linked. 
We then use the element-by-element analysis in order to discuss 
our results compared to the current literature on the halo-disk 
transition regions in section 4. In section 5 we use the abundance 
ratio trends to propose a new powerful chemical labeling technique 
that uses only chemistry to efficiently separate the various Galactic 
components. Finally, we summarize our findings in section 6. 

2 DATA AND SUBSAMPLES 
2.1 Data: The APOGEE Survey 

To study the interface of the Galactic disks and halo, we employed 
the first two years of data from the SDSS-III/APOGEE survey (de¬ 
scribed in Eisenstein et al. 2011). APOGEE used H band spec¬ 
troscopy to determine the stellar parameters, chemical abundances. 


and radial velocities for a large sample of stars. The most cur¬ 
rent data release (SDSS DR12) contains nearly 100 000 stars with 
stellar parameters, chemical abundance determinations and spectra 
(Holtzman et al. 2015). The typical uncertainties in the parame¬ 
ters are ±100 K, ±0.1 dex, ±0.04 dex, and ±0.04 dex in TeS, 
log g, [Fe/H], and [a/Fe], respectively. While there may be some 
systematics biases (potentially up to -1-0.20 dex, see Masseron & 
Gilmore 2015), the precision of APOGEE chemical abundances are 
expected to be relatively high (with typical internal uncertainties of 
less than 0.10 dex). To minimize the potential uncertainties in the 
accuracy in the abundance scales, in this paper we have done a rel¬ 
ative abundance analysis of the various Galactic components. 

We selected a sample of giant stars, beginning with the full 
APOGEE sample and using the following cuts: 

• The signal-to-noise ratio (SNR) had to be larger than 100. We 
used this cut to ensure that the selected stars have quality estimates 
of the stellar parameters and [a/Fe]. 

• Teff > 4000 K. This cut was chosen because as noted by 
Holtzman et al. (2015), stars with Tgff < 4000 K will likely have 
larger uncertainties in stellar parameters and chemical abundance. 

• The ASPCAR metallicity, and [a/Fe] flags must be set to 0. 
This cut ensured that there were no major flagged issues (e.g. low 
SNR, poor synthetic spectral fit, stellar parameters near grid bound¬ 
aries, etc.). 

• 1.0 < log g < 3.5 dex. This cut was used to deselect any 
dwarf stars for which the stellar parameters and chemical abun¬ 
dances would likely have large uncertainties (e.g. Nidever et al. 
2014; Holtzman et al. 2015). 

• The stars had to be located at 6 < 60 degrees. This was needed 
to separate disk-like from non-disk like kinematics in /-Galactic 
rest frame radial velocity (GRV) space as we discuss in section 2.2. 
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The initial sample contained about 69000 stars. The sample 
was then restricted to the metallicity regions -1.20 < [Fe/H] < 
-0.55 dex in which encompasses the transition between the Galac¬ 
tic disks and halo. This reduced the sample to about 3200 stars. We 
note that we used the uncalibrated [Fe/H] for metallicity rather than 
the calibrated global metallicity, [M/H], for the purposes of com¬ 
paring [X/Fe] abundance ratios to literature. Holtzman et al. (2015) 
noted that the APOGEE uncalibrated [Ee/H] tends to be overes¬ 
timated at low metallicities. However, applying the calibration to 
[Ee/H] would cause the [X/Ee] ratios to be inflated compared to 
literature values. In this paper we focus on a relative abundance 
analysis within subgroups, and as such the use of [Fe/H] as op¬ 
posed to [M/H] only changes the absolute abundance scale but does 
not change the relative analysis leaving our conclusions robust. See 
section 4.3 for more details. In the following section we use the 
chemistry and kinematics to decompose the sample into the vari¬ 
ous canonical Galactic components (similar to the scheme outlined 
in Figure 1). 

2.2 Chemokinematic decomposition of canonical Galactic 
components 

Since we are primarily interested in the disks-halo interface, our fo¬ 
cus in this paper is on the intermediate metallicity (-1.20 < [Fe/H] 
< -0.55 dex) region. In Figure 3(a) we plot the standard [a/Fe]- 
metallicity for the intermediate metallicity stars. With our initial 
sample we recovered the two disk sequences discussed in section 

1.2 and seen in other APOGEE studies (e.g Nidever et al. 2014; 
Masseron & Gilmore 2015). In the low-metallicity regime (-1.20 < 
[Ee/H] < -0.70 dex), we found the two sequences: an a-poor se¬ 
quence (left magenta selection box) and an a-rich sequence (the 
left orange selection box) similar to Nissen & Schuster (2010). We 
determined the boundary between these two sequences by identi¬ 
fying the trough in the [a/Fe] distribution in three metallicity bins 
and performing a linear fit in the [a/Fe]-metallicity plane of these 
three metallicity bins similar to Recio-Blanco et al. (2014). In Fig¬ 
ure 2, we show the [a/Fe] distribution in each of the metallicity 
bins used to determine the separation between the a-poor and a- 
rich sequences. The bins were chosen to ensure at least 50 stars 
were in each bin. We note that the bimodality is most clear when 
all halo stars (i.e. those stars assigned only to canonical and ac¬ 
creted halos not including stars from the disks and undetermined 
groups) are used without binning in metallicity space. 

With the lack distances and proper motions we do not employ 
the Toomre diagram as a means assign the stars to a component be¬ 
cause we cannot estimate their full 3D velocity. However, we made 
use the RVs, which have typical uncertainties on the order of 0.1 
km s“^. We used a method that combines the I and GRV in order 
to disentangle the halo from disk stars similar to Sheffield et al. 
(2012). We computed the GRV by correcting the radial velocity 
for solar motion (see Equation 1 of Hawkins et al. 2015). Eor this, 
we assumed vq = [Uq, Vq,Wq\ = [14.0,12.24, 7.25] km s“^, 
Vlsr = 220 km s“^ (Schonrich 2012). In Eigure 3(b) we show 
the GRV/cos(b) as a function of I for the full APOGEE sample 
that pass our quality cuts (black contours). As noted by Majewski 
et al. (2012) and Sheffield et al. (2012), dividing by the GRV by 
the cos(6) accentuates the difference between the disk and non-disk 
components provided the data is restricted to 161 < 60 degrees al¬ 
lowing us to separate out halo stars more efficiently. This cut, while 
potentially eliminating some halo stars, still leaves enough to study 
the chemical abundance trends. We have also tested several 161 cuts 
between 40 - 60 deg with little to no affect on the chemical abun- 



[a/pq 


Figure 2. The distribution of [a/Fe] in three metallicity bins, with bin cen¬ 
ters at [Fe/H] = —0.96, —0.85, —0.76 dex respectively. In each metallicity 
bin there appears to be a bimodal distribution in [a/Fe] the trough between 
the a-poor and a-rich components are marked by a black arrow. 

dance trends and overall conclusions of this work. As 161 approaches 
larger values, the cos(6) becomes too small relative to the GRV to 
be useful. 

In /-GRV space (Figure 3(b)), the Galactic disk makes a sinu¬ 
soidal pattern with an amplitude approximately equal to the rota¬ 
tional velocity of the Galaxy and dispersion correlated to the ve¬ 
locity dispersion of the disk. In the same space the Galactic halo 
has a more-or-less random distribution (e.g. Johnston et al. 2012). 
To characterize the disk in this space, we selected a subsample of 
APOGEE stars with metallicities between -0.50 and 0.00 dex (i.e. 
the metallicity range where the disk dominates) and binned in I 
space (20 degrees bins). The dotted blue lines in the top panel of 
Figure 3(b) are equal to 2.5 times the GRV dispersion in each I bin 
which represents stars with disk-like kinematics. Stars that have an 
absolute GRV/cos(6) that is at least 2.5a or larger than the Galac¬ 
tic disk within a given I bin are considered to have non-disk like 
kinematics and are assigned to the halo. We note that there are a 
handful of stars which are assigned to the Galactic halo which ap¬ 
pear to sit below GRV dispersion cut in Figure 3(b). These stars 
are on the I bin edges and assigning them to the halo or disk does 
not impact the results given the large sample size. While there still 
may be some disk contamination in the halo sample, this is likely 
to be small. Similarly, there may be halo stars with disk-like kine¬ 
matics that this study also rejected. Full 3D space motions could 
aid in recovering these lost halo stars while further rejecting thick 
disk stars. We varied the GRV dispersion cut from 1 - 4 a to asses 
the impact of the kinematic classification. We discuss the impact of 
the kinematic selection in section 4, however we note here it does 
not affect our conclusions. 

As we are interested in studying the bulk properties of the 
thin-thick disk-halo transition, we must use the chemodynamical 
information to define our sample into the canonical Galactic com¬ 
ponents discussed in sections 1.1 and 1.2. Our selection is inspired 
by the literature described in section 1.1, however it is important 
to keep in mind that our definitions may not be the same as other 
studies thus comparisons should be done with care. We define the 
Galactic components in our sample in the following way: 
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Figure 3. (a) The [a/Fe] as a function of metallicity for the intermediate metallicity APOGEE sample. We separated chemical the various Galactic components 
based on where they sit in [a/Fe]-metallicity diagram. A cut is placed at [Fe/H] = -0.70 dex as this where the thin disk is thought to end (e.g. Bensby, Feltzing 
& Oey 2014). The error bar represents the median uncertainty in [o/Fe] and metallicity. The black points are the stars that were deselected because it is difficult 
to classify those stai's on the basis of their chemistry, (b) Upper panel: The GRV as a function of I for the full APOGEE sample that pass our quality cuts 
(contours). The dotted blue lines represents 2.5 x disk GRV dispersion in each I bin. It is clear from this diagram that the disk signature is a sinusoidal pattern 
in /-GRV space. Therefore we select stars that have absolute GRVs larger than dotted blue line as our non-disk subsamples (i.e. these stars have, on average 
less than a 1.24 per cent chance of being kinematically a part of the disk), which are marked as black squares and red triangle for the a-poor and a-rich, 
respectively. While the orange stai's and cyan diamonds ai‘e the a-rich and a-poor stars that have disk-like kinematics (i.e. they are inside of the blue dotted 
line) and thus are consistent with disk-like kinematics. Stars with metallicities larger than -0.70 dex that are a-poor are considered canonical thin disk stai's 
(open magenta hexagons) Lower Panel: the /, b Galactic coordinates of the full APOGEE sample that pass our quality cuts (contours), following the same 
symbols as the top panel. 
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• Canonical Thin Disk (open magenta hexagon) - These are 
stars (shown as the right magenta selection box in Figure 3(a)) that 
have low [a/Fe] ratios ([a/Fe] < +0.10 dex, e.g. Adibekyan et al. 
2012; Haywood et al. 2013; Recio-Blanco et al. 2014) and have 
metallicities > -0.70 dex (e.g. Bensby, Feltzing & Oey 2014). The 
cut in [a/Fe] was designed to conservatively select out the Galactic 
thin disk. Kinematically these stars have spatial-kinematic coher¬ 
ence (falling inside of the blue dotted line in Figure 3(b)). We note 
that a result of cutting our sample at [Fe/H] < -0.55 dex and the 
negative radial metallicity gradient found in the disk (e.g. Cheng 
et al. 2012), the thin disk stars in our sample are preferentially 
found in the outer Galaxy. 

• Canonical Thick Disk (orange stars) - These stars (orange 
selection boxes in Figure 3(a)) have high [a/Fe] ratios ([a/Fe] 
larger than +0.22 dex at [Fe/H] > -0.70 dex, e.g. Adibekyan et al. 
2012; Haywood et al. 2013; Recio-Blanco et al. 2014). Kinemati¬ 
cally these stars have spatial-kinematic coherence (falling inside of 
the blue dotted line in Figure 3(b)). 

• Canonical Halo (black squares) - These stars are defined as 
a-rich, metal-poor stars (left orange selection box in Figure 3(a)). 
The [a/Fe] cut was determined by a linear fit to the trough in the 
[a/Fe] distribution in three metallicity bins below -0.70 dex, where 
the halo and thick disk will dominate. Kinematically they were se¬ 
lected to have high velocities (e.g. Nissen & Schuster 2010; Schus¬ 
ter et al. 2012; Hawkins et al. 2015) and not have significant spatial- 
kinematic coherence (falling outside of the blue dotted line in Fig¬ 
ure 3(b)). Additionally, there may be some contamination of this 
component from the accreted halo caused by the ‘hard’ cut used to 
separate the two. However, less than a fifth of the canonical halo 
stars fall within \-a of the [a/Fe] boundary. 

• Accreted Halo (red triangles)- These stars are a-poor and 
metal-poor (left magenta selection box in Figure 3(a)). They were 
selected kinematically to have high velocities (e.g. Nissen & Schus¬ 
ter 2010; Schuster et al. 2012) and not have significant have spatial- 
kinematic coherence (falling outside of the blue dotted line in Fig¬ 
ure 3(b)). Our choice to use a ‘hard’ boundary to separate the a- 
poor accreted halo and the a-rich canonical halo (as done in other 
literature, e.g. Nissen & Schuster 2010, 2011; Schuster et al. 2012; 
Bensby, Feltzing & Oey 2014) may result in some contamination 
of the accreted halo subsample from the a-rich canonical halo (or 
thick disk in some cases) or vice-versa. Accreted halo stars within 
~ 0.02 dex of the boundary, the typical \-a uncertainties in [a/Fe], 
could potentially be true canonical halo stars. This only accounts 
for less than a quarter of the total accreted halo sample. Thus the 
canonical halo is not likely to be a large contaminating source in 
the accreted halo. 

• Undetermined (cyan diamonds)- These stars are a-poor 
(right magenta selection box in Figure 3(a)) with metallicities less 
than -0.70 dex but may have spatial-kinematic coherence (falling 
inside the blue dotted line in Figure 3(b)). In this way these stars 
have the chemistry of what would be the accreted halo but kinemat¬ 
ics that cannot rule out disk-like motion. The lowest metallicity thin 
disk stars have been found to have [Fe/H] = -0.70 dex (e.g. Bensby, 
Feltzing & Oey 2014), thus these undetermined stars could be a 
part of the accreted halo or misidentified thick disk stars with lower 
[a/Fe]. Additionally the undetermined stars within ~ 0.02 dex of 
the boundary separating a-rich and a-poor stars may be canonical 
halo stars however, the lower the [a/Fe] of the undetermined stars 
the less likely that this is the case. 

The halo is thought random in Z-GRV space (e.g. Johnston 
et al. 2012) and thus with a lack of full 3D space motions there is 


a degeneracy that could result in some contamination of halo stars 
in all canonical disk subgroups. However, even though the halo is 
thought to extend to metallicities as large as -0.50 dex (e.g. Nissen 
& Schuster 2011), the Galactic disk populations will dominate at 
metallicities above -0.80 dex (e.g. Schlesinger et al. 2012) particu¬ 
larly if they follow the sinusoidal disk-like pattern in Figure 3(b). A 
rough estimate of the upper limit of the contamination of the disk 
components on the halo components as a result of using the /-GRV 
diagram to decompose the Galaxy can be determined using the pro¬ 
cedure outlined in section 5 of Sheffield et al. (2012). Assuming 
that all 3200 stars in our sample are disk stars, and that the stars 
in a given I bin are distributed in a Gaussian way around the un¬ 
derlying sinusoidal ‘disk-like’ feature, ~ 40 of those stars (i.e. 1.2 
per cent or 2.5 a) would be identified as ‘halo’ stars. However, we 
found nearly 160 potential halo stars of which we estimate at most 
40 of these could be true disk members. This would yield a con¬ 
tamination rate of 25 per cent. Additionally, we used the Galaxia 
Galactic model (Sharma et al. 2011) to further study the expected 
contamination rate. The models have been set to have a metallicity 
distribution for each component described in section 3.1 of Kor- 
dopatis et al. (2013). We used the same magnitude and color cuts 
of APOGEE. We applied the same selection routine as above. Re¬ 
sults indicate that the contamination rate of the (thick) disk on the 
Galactic halo samples can range between 20-60 percent depend¬ 
ing on the maximum distance surveyed. Larger sampling distance 
yield lower contamination rates. We also note that adopting a 3 -(j 
cut, instead of 2.5-a, reduces this contamination to in most cases 
below 15 per cent. As such, we varied the GRV dispersion cut to as 
high as 4-a (changing the amount of contamination in the disk and 
halo subgroups) and found that the chemical abundance trend and 
the overall conclusions are not significantly affected. 

Our final sample contains a total of approximately 150 unde¬ 
termined stars (cyan diamonds) in the low-a, low-metallicity re¬ 
gion with disk-like kinematics that could be thin disk stars or ac¬ 
creted stars. There are also approximately 160 halo-like stars split 
in almost equal parts between the accreted halo (red triangles) and 
the canonical halo components (black squares). The accreted halo 
seem to be distributed over all I, b space that has been sampled. 
Unsurprisingly there are more halo stars at high Galactic latitudes. 
Accurate distance and proper motions could help compare the total 
space motions and orbital parameters of the various components, 
which is beyond the scope of this paper. There are approximately 
370 canonical thin disk stars (open magenta hexagons), and 1400 
canonical thick disk stars (orange stars). 


3 ANALYSIS OF CHEMICAL ABUNDANCES OF THE 
CANONICAL GALACTIC COMPONENTS 

In this section we discuss the individual abundance ratios and 
trends as a function of metallicity of the five identified subsamples 
(i.e. canonical thin disk, canonical thick disk, halo and accreted 
halo and undetermined) defined in section 2.2. Although one may 
be concerned by the systematic errors of the absolute abundances 
in the APOGEE data, our study mainly relies on relative compar¬ 
ison between the different populations. Therefore, our discussion 
and conclusions should not be affected by these systematic errors; 
in our case only the relative abundance precision matters. 
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3.1 The a-elements: O, Mg, Si, S, Ca, and Ti 

We study the a-elements as a way to infer the properties of the 
size and star formation rate of the cloud that formed the different 
stellar populations. The a-elements are those elements that are, by 
definition, made by capturing a-particles. This implies that those 
elements have an atomic mass number that is a multiple of 4, 
thus Mg,^® Si,®^ S,'*° Ca, and '^"‘Ti. While other isotopes 

of those elements exist, they cannot currently be distinguished in 
the APOGEE dataset, as such we will assume that the listed iso¬ 
topes are the dominant contributor to the elemental abundances of 
the APOGEE data. 

Those elements are made in the cores of stars by a-capture 
during post-main sequence burning and dispersed in the interstellar 
medium via type II supemovae (SNII), while relatively less is ex¬ 
pected to be produced by type la supernovae (SNIa). Thus the ratio 
of [a/Ee] is thought to be sensitive to the environment enrichment 
history of the gas from which the star formed. Unlike the a-poor 
stars, the a-rich stars are thought to be formed in regions where 
SNII were numerous, thus in a significantly higher star formation 
than the a-poor stars (e.g. Gilmore & Wyse 1998). 

In Eigure 3(a) we see that indeed APOGEE data shows evi¬ 
dence of an a-poor and a-rich sequence at both the intermediate 
and low metallicity regimes, indicating a distinct formation his¬ 
tory for each sequence (Nidever et al. 2011). In Figure 4 we plot 
the abundance of each a-element relative to iron as a function of 
metallicity in a different panel, ordered by atomic number. The dif¬ 
ferent symbols represent the populations defined in section 2.2. The 
abrupt change of symbols at [Fe/H] = -0.70 dex is due to the defi¬ 
nition of the populations only. 

Unlike several Galactic studies in the literature (e.g. Reddy 
et al. 2003; Venn et al. 2004; Haywood et al. 2013; Bensby, Feltzing 
& Oey 2014), our sample contains the a-element sulphur. Not sur¬ 
prisingly, it behaves in a similar way to the other a-elements (e.g. 
Nissen et al. 2007; Spite et al. 2011; Jdnsson et al. 2011; Takada- 
Hidai & Takeda 2012; Matrozis, Ryde & Dupree 2013; Skuladottir 
et al. 2015) yet not used in the standard [a/Fe] definition. Since we 
used this the global [a/Fe] criteria to define our populations, we ob¬ 
viously find that our populations separate in each of the individual 
a-elements. 

From Figure 4 we also note that in the -0.90 < [Fe/H] < - 
0.70 dex metallicity regime, O and Mg drive the largest distinction 
between the a-rich and a-poor sequences. The rest of the a el¬ 
ement abundances appear to be very mixed, making it difficult to 
disentangle a clear a-poor sequence from an a-rich one. This could 
be a result of O and Mg SNII yields being mass-dependent (e.g. 
Nomoto, Kobayashi & Tominaga 2013). 


3.2 The Fe-peak elements: Mn, Ni 

In contrast to the a-elements, Fe-peak elements, such as Mn and 
Ni, are thought to be formed primarily via SNIa (e.g. Iwamoto et al. 
1999). Therefore they are expected to track Fe in such a way that 
[Mn/Fe] has been shown to be sub solar in low metallicity stars 
and increase as with [Fe/H] (e.g. Adibekyan et al. 2012; Battistini 
& Bensby 2015). We plot in Figure 5 the abundances of the only 
Fe-peak elements in APOGEE, namely Mn and Ni, for the inter¬ 
mediate metallicity stars as a function of [Fe/H]. As in Figure 4, 
the different symbols represent our definition of populations (see 
Sect. 2.2). 

It is interesting to note that Mn has a different behaviour for 
thin and thick disk stars (orange and magenta colour, respectively) 


at comparable metallicities. While thin disk stars have Mn en¬ 
hanced respect to iron, thick disk stars have it at the same level 
than iron. Since Mn is produced at a higher fraction compared to 
Fe during SNIa (e.g. Gratton 1989; Iwamoto et al. 1999; Kobayashi 
et al. 2006) one may expect that at a given metallicity, a-poor 
stars - which have been polluted by more SNIa - will have higher 
[Mn/Fe] ratios compared to their a-rich counterparts. This offers 
an explanation as to why the canonical thin disk stars seem to have 
overall higher [Mn/Fe] ratios compared to the canonical thick disk. 
Kobayashi & Nakasato (2011) demonstrate from models that O and 
Mn should be anti correlated in the Galaxy, because O is produced 
mainly by SNII while Mn is produced predominantly by SNIa. We 
observe in Figure 4 that the canonical thin disk has lower [0/Fe] 
than its thick disk counterpart. Thus, given the results of Kobayashi 
& Nakasato (201 1), one would expect the thin disk to have a higher 
[Mn/Fe] ratio than the thick disk. 

The same authors also predict Ni to behave similarly to Mn. 
However, in Figure 5 the [Ni/Fe] remains constant between thin 
and thick disk stars. More puzzling is the fact that [Ni/Fe] behaves 
the opposite at low metallicities, in which the Ni abundances of the 
a-poor accreted population are lower than the those of the a-rich 
halo at same metallicities. This agrees with the results of Nissen 
& Schuster (2010), but there is currently no explanation for this 
unexpected tendency. 

We note in Figure 5 the very tight dispersion of stars within 
each population (i.e. stars which have the same symbols). This con¬ 
firms that despite APOGEE data requiring some further calibra¬ 
tions of the absolute values of the elemental abundances, they pro¬ 
vide good precision for Mn and Ni, demonstrating that this dataset 
is excellent for chemical tagging purposes. Therefore they are very 
good candidates to allow accurate distinctions between our popu¬ 
lations, besides the fact that the comparison of the absolute abun¬ 
dances of [Mn/Fe] and [Ni/Fe] with external studies might be non 
trivial. 

3.3 The Odd-Z, Even-Z, and Light Elements 

C+N 

APOGEE stars are giants, and as such their surface C and N abun¬ 
dances have been partially affected during their evolution as a result 
of dredge up processes (e.g. Iben 1965). However, because a large 
majority of the stars of this survey are expected to be low-mass stars 
(« IMq), where CN cycling processes occurred, the initial C-l-N 
ratios are conserved throughout the evolution of those stars (see 
Masseron & Gilmore 2015). Therefore, we can use C-l-N to discuss 
and compare the evolution of C in our defined populations. In Fig¬ 
ure 6 we display in top panel [C-l-N/Fe] as a function of metallicity 
for our defined populations with the same symbols as previously. 
Carbon is mostly made by He burning and its main contributors 
into the enrichment of the Galaxy are SNII at very low metallic¬ 
ity as well as AGB stars at around -1.50 dex in metallicity. Hence, 
C-l-N is expected to globally increase with metallicity, but decrease 
as soon as SNIa kicks in because it does not produce C. This is 
why in Figure 6, we can observe some enhancement of [C-l-N/Fe] 
for the lowest metallicity stars, but then there is an overall decrease, 
in particular in the thin disk which is the population most enriched 
by SNIa. Because of this unique interplay among other elements 
between different progenitors, and thanks to the good precision of 
C-l-N in the APOGEE data, the [C-l-N/Ee] ratios offer an opportunity 
to disentangle Galactic populations. 
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Figure 4. The [0/Fe], [Mg/Fe], [Si/Fe], [S/Fe], [Ca/Fe], [Ti/Fe] from top to bottom respectively as a function of metallicity for the accreted halo (red 
triangles), canonical halo (black squares), canonical thick disk (orange stars) and undetermined (cyan diamonds), and canonical thin disk (magenta open 
hexagons) subgroups. It is clear that the distinction between the a-rich and a-poor sequence at low metallicity is likely driven by the [Mg/Fe] and [0/Fe] 
abundance ratios. 
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Figure 5. The [Mn/Fe], [Ni/Fe] from top to bottom, respectively as a function of metallicity for the five groups with the same symbols as Figure 4. 


Odd-Z elements: Na, Al 

The second and third panels of Figure 6 show the trends of [Na/Fe] 
and [Al/Fe] with metallicity, with the symbols representing the 
populations described in Sect. 2.2. While Na abundances are rel¬ 
atively scattered at low metallicities, Al tends to separate. It is ex¬ 
pected the Na and Al are produced in SNII. However, according 
to Kobayashi et al. (2006), the production of those elements are 
strongly dependent on the initial C and N in the gas cloud that forms 
the stars. Therefore, it is expected that Na and Al are correlated with 
C-fN. Besides the SNII production site, Na and Al are also expected 
to be partially produced by AGB stars in the Galaxy as shown by 
Nomoto, Kobayashi & Tominaga (2013) at the metallicities we fo¬ 
cus on in this paper. This explains why in Figure 6 [Al/Fe] is rather 
enhanced in the lowest metallicity part of the thick disk population, 
because [C-l-N/Fe] is also enhanced. Moreover since SNIa do not 
produce Al and Na as efficiently as Fe, the [Al/Fe] and [Na/Fe] tend 
to decrease towards higher metallicities. Nissen & Schuster (2010) 
demonstrate the effective ability of Na to characterize the a-poor 
populations. Unfortunately at these metallicities APOGEE has dif¬ 
ficulty measuring Na (see Figure 14 of Holtzman et al. 2015). Typ¬ 
ical [Na/Fe] abundance uncertainties in this metallicity range ap¬ 
proach ~ ± 0.20 dex. In contrast, the precision of Al abundances 
in the APOGEE data is effectively very high, and thus offers an 
alternative to Na. 

KandV 

The last bottom panels of Eigure 6 display the trends of potassium 
and vanadium as a function of metallicity for our different defined 
populations with different symbols. The abundance V shows 
a very large dispersion and remains rather constant among all 
populations, thus tracking Fe. This is similar to what is seen in the 
literature for V (e.g. Bodaghee et al. 2003; Battistini & Bensby 
2015). Concerning K, the abundances of this element seem to 
decrease towards higher metallicity consistent with literature (e.g. 
Shimansky et al. 2003). One can also notice a decrease in some 
of the accreted population stars, but this trend is rather shallow. 
Although we note that the nucleosynthesis channel is rather 
misunderstood for K and V, and the supernovae yields lead to 


largely underestimated values compared to their observed Galactic 
abundances (e.g. Nomoto, Kobayashi & Tominaga 2013). As such, 
these elements do not give strong constrain to distinguish our 
defined populations. 

With the element-by-element discussion above, in the next 
section we focus on using the abundance ratio trends to inform 
our understanding of the transition between the Galactic disks and 
halo. 


4 IMPLICATIONS OF CHEMICAL ABUNDANCE 
TRENDS 

In the previous section, we have demonstrated that among a- 
elements, O and Mg bear distinction of the star formation activity, 
particularly at low metallicities. It also indicates that Mn and Ni are 
potential alternatives to disentangle the signatures of SNIa. Regard¬ 
ing the accreted population versus thin disk, Al and C-l-N represents 
a good diagnostic as they contain signatures of AGB nucleosynthe¬ 
sis. The observations of those elemental abundances provide insight 
about their origin, but this requires some further modeling taking 
into account stellar yields as well as a thoughtful discussion of the 
observational bias that is beyond the scope of this paper. Neverthe¬ 
less, we will rather use their relative abundance ratios in this section 
to discuss the connection and the transition phases of the different 
canonical components of the Milky Way. 

4.1 Thick disk-halo transition 

One of the primary goals of this work is to use the chemistry of a 
large sample of stars to study the disk-halo transition in the metal¬ 
licity region where they are thought to overlap (-1.20 < [Fe/H] 
< -0.70 dex) based on literature (see section 1.2 for more de¬ 
tails). In each element studied, the canonical halo (black squares) 
and canonical thick disk (orange stars) components are completely 
indistinguishable (see Figures 4, 5, 6). This result was tested us¬ 
ing different GRV dispersion cuts between I -4a (i.e. varying the 
amount of contamination between the two groups) and in all cases 
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Figure 6. The [C+N/Fe], [Na/Fe], [Al/Fe], [K/Fe], and [V/Fe] as a function of metallicity for the five subgroups are shown with same symbols as Figure 4. 


the canonical thick disk and halo overlap in every element. Another 
illustration of this can he seen in Figure 7 for the key elements we 
have previously highlighted. In each panel of that figure the canon¬ 
ical thick disk (shown in orange) and the canonical halo (shown in 
black) not only overlap but have comparable dispersions suggest¬ 
ing that the two formed in a similar high star formation environ¬ 
ment with well-mixed turbulent gas making them indistinguishable 
in chemistry leaving only the kinematics and spatial distributions 
to separate the two. 

Ishigaki, Chiba & Aoki (2012) found, using a sample of 97 
stars, evidence that [Mg/Fe], [Si/Fe] in the thick disk are higher, on 
average, than halo stars. Their follow-up study (Ishigaki, Aoki & 


Chiba 2013) found that [Na/Fe] and [Ni/Fe] were lower in the halo 
compared to the thick disk. It is important to note that we select our 
halo samples in different ways making it difficult to compare these 
two studies. The halo sample in Ishigaki, Chiba & Aoki (2012) 
was selected kinematically on the basis of the Toomre diagram and 
is split into an ‘inner halo’ and ‘outer halo’ component. It is likely 
that both of their inner and outer halo components contain stars that 
we would classify as a-rich and a-poor. In our study, we selected 
the halo sample kinematics (using the Z-GRV diagram) and found 
a bimodal distribution in [a/Fe]. As a result of this we called the 
a-rich components the ‘canonical halo’ and the a-poor component 
the ‘accreted halo’ (based on studies such as Nissen & Schuster 
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2010). Since both halo samples in Ishigaki, Chiba & Aoki (2012) 
likely contain a-rich and a-poor stars, the mean [ct/Fe] could be 
driven to lower values which may explain why the [Mg/Fe] and 
[Si/Fe] may have been found to be on average lower in the halo(s) 
than the thick disk. 

The only chemical distinction between thick disk and canon¬ 
ical halo is solely in the mean and dispersion of the metallicity 
distribution function for the two components. This result could be 
biased by the fact that we choose the canonical halo to be a-rich, 
similar to the thick disk. It is interesting to point out that while we 
selected the canonical halo and thick disk to be a-rich, they are, 
in fact, similar in all other elements including elements that have 
different production channels to [a/Fe] (e.g. odd, even, and light 
elements) that we do not use for selection. If we had not initially 
separated the halo sample into a-rich and a-poor subsamples, we 
would conclude that in key elemental abundance spaces the halo is 
bimodal (see Figure 7) and the a-rich component of that bimodal 
distribution is chemically similar to the thick disk despite having 
a different metallicity distribution function while the a-poor com¬ 
ponent is chemically different. It is also interesting to note that the 
canonical halo has some evidence of spatial-kinematic coherence 
(black squares in Figure 3(b)), indicating it may have net rotation 
in line with other studies (e.g. Carollo et al. 2007, 2010). This may 
suggest that the thick disk-halo transition is smooth kinematically 
and would be worth examining with full 3D space motions to con¬ 
firm. As the thick disk is thought to form from well-mixed tur¬ 
bulent gas (e.g. Haywood et al. 2013), we interpret this smooth 
chemical transition between the canonical thick disk and halo com¬ 
ponents as evidence that the gas of the inner in situ Galactic halo 
may have been the precursor to the thick disk. We note that subject 
to careful modeling of selection function, most of the inner halo 
(i.e. the a-rich stars) chemically resembles a single homogeneous 
self-enriching population, in which the major transition between 
the halo and thick disk is pressure support to angular momentum 
support in kinematic space. 

4.2 The accreted halo 

In section 2.2 we showed that there is a sequence of a-poor metal- 
poor halo giant stars in APOGEE (which we designated as ‘ac¬ 
creted halo’) that resembles the a-poor sequence of Nissen & 
Schuster (2010), Ramirez, Melendez & Chaname (2012), Bensby, 
Feltzing & Oey (2014), Hawkins et al. (2014), and others. What 
makes this result interesting is that APOGEE has a very different 
selection function and surveys a larger volume than past studies 
(e.g. Nissen & Schuster 2010) which primarily rely on chemical 
analysis of relatively local dwarf stars. The discovery of the a-poor 
sequence in a large sample of giant stars from APOGEE with a dif¬ 
ferent selection function and larger survey volume strengthens the 
argument that the a-poor halo sequence is real. 

We used the APOGEE survey to study the chemistry of these 
a-poor stars (red triangles) in section 3 by comparing their abun¬ 
dances in 14 elements, including a (Figure 4), Fe-peak elements 
(Figure 5), odd, even, and light elements (Figure 6) to the canonical 
halo sequence (black squares). We found that the accreted halo is 
primarily distinguishable in the a-elements (e.g. Mg, O, and S). To 
further illustrate the chemical distinction between the accreted halo 
and canonical halo sequences, in Figure 7 we plotted the chemical 
distribution of six key elements from all three elemental groups for 
the canonical halo (black line), accreted halo (red line) components. 

In all chemical elemental abundance ratios shown in Figure 
7 except [Mn/Fe] the accreted halo is systematically under abun¬ 


dant compared to the canonical halo. The lower [C-l-N/Fe] of the 
accreted halo indicates that the initial chemical composition of the 
gas that formed the a-poor stars were different from the a-rich pro¬ 
viding further evidence that the a-poor stars are likely extragalac- 
tic in origin (e.g. Nissen et al. 2014). This is the first observation 
that [Al/Fe] is under abundant in the accreted halo population. This 
result is not surprising as Al is thought to be produced in SNII 
and correlated strongly to C-l-N. Therefore, if [C-l-N/Fe] is under 
abundant (see section 3.3 and Nissen et al. 2014) one would expect 
[Al/Fe] to also be under abundant. Furthermore, the dispersion of 
the accreted halo in all elemental ratios except Mn is larger than the 
canonical halo component. This may suggest several extragalactic 
systems being responsible for the production of the accreted halo 
component in contrast to the canonical halo, which is thought to 
form from one homogenous gas cloud. However, we remind the 
reader that the dispersion in each element is sensitive to the bound¬ 
ary that separates the canonical and accreted halo components. As 
such, we tested various boundaries between the canonical and ac¬ 
creted halo components to study how this would affect the disper¬ 
sions and the separation of the two components. We shifted the 
line that separates the a-rich canonical and a-poor ‘accreted’ ha¬ 
los by at least ±0.05 dex (twice the typical internal uncertainty in 
[a/Fe]). Lowering the boundary causes stars currently classified as 
accreted with moderate [a/Fe] to be reclassified as canonical halo. 
The overall result is a decrease in the dispersion of accreted halo in 
most elemental spaces and an increase in the dispersion in the el¬ 
emental distributions for the canonical halo. Raising the boundary 
has the opposite effect. In both cases, the accreted halo is chemi¬ 
cally distinct from the canonical halo in the elements listed above. 
If we were to not separate the halo into a ‘canonical’ and ‘accreted’ 
component at the onset, we would find a bimodal distribution in O, 
Mg, Al, C-l-N and Ni in the full halo sample consistent with two 
chemically distinct components. 

As noted by Nissen & Schuster (2011) the [Mn/Fe] ratio is 
comparable between the a-poor and a-rich halo sequences. This 
is likely a result of the metallicity dependence on the SNIa yields. 
We also found that [K/Fe] and [V/Fe] ratios are comparable both 
in mean and dispersion between the accreted and canonical ha¬ 
los indicating that these elements do not give strong constraints 
to distinguish the canonical and accreted halo components. We 
recall that the [V/Fe] ratios have significantly larger scatter than 
expected from other literature (e.g. Fulbright 2000; Battistini & 
Bensby 2015) which may indicate this element is less reliable. 

Additionally, increasing the kinematic GRV dispersion cut to 
4a reduces the contamination of the disk population in the halo 
stars however decreases the sample size significantly while reduc¬ 
ing the kinematic GRV dispersion cut to Icr increases the con¬ 
tamination from the a-poor disk population thereby broadening 
the chemical abundance ratio dispersions of the a-poor ‘halo’ se¬ 
quence. In all cases in which we tested various GRV dispersion 
cuts, the canonical and accreted halo remain chemically distinct 
in the [C-l-N/Fe], [Al/Fe], [Ni/Fe] and the a-elements (particularly 
[O/Fe], [Mg/Fe], [S/Fe]). 

Simulations of the Galactic halo (e.g. Bullock & Johnston 
2005; Zolotov et al. 2009; Cooper et al. 2010) have indicated that 
the halo may have been constructed primarily through a combi¬ 
nation of dissipative collapse and accretion events of massive sys¬ 
tems. The likely accreted a-poor stars represent accretion from sys¬ 
tems onto the halo and their [a/Fe] and chemical abundance distri¬ 
butions may be able to give insight to the mass of the systems that 
accreted. Furthermore, a careful study of the relative fraction of 
the accreted a-poor halo stars and the in situ a-rich halo stars will 
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Figure 7. The distribution of [0/Fe], [Mg/Fe], [Al/Fe], [C+N/Fe], [Ni/Fe], [Mn/Fe] abundance ratios for the accreted halo (red), canonical halo (black) and 
canonical thick disk sequence (orange). We find that in most elements (except Mn, Na) are systematically under abundant in the accreted halo sequence relative 
to the canonical halo sequence. The abundance dispersion in the accreted halo sequence is comparable or in most cases larger than the canonical halo sequence. 
Further Mn shows no chemical difference between the two sequences. In all chemical abundance ratios shown (and those not shown as well) the canonical 
thick disk and halo components overlap significantly. 


provide unique constraints on the relative importance of accretion 
events in the assembly of the Galactic halo. With the upcoming sur¬ 
veys it will be possible to build large samples of a-poor halo stars 
that will make these tests possible. 

4.3 The undetermined group: thin disk at low metaUicitles? 

The thin disk is thought to be rather metal-rich with [Fe/H] > 
—0.70 dex (e.g. Reddy et al. 2003; Fuhrmann 2004). However, 
there has been some debate as to whether the thin disk actually ex¬ 
tends to metallicities as low as [Fe/H] = -1.0 dex (e.g. Mishenina 
et al. 2004). Most recently Bensby, Feltzing & Oey (2014) used a 
sample of 714 F and G dwarf stars and found that the lower metal- 
licity for the thin disk was [Fe/H] ~ -0.70 dex. Their finding in 
combination with the extent of the accreted halo (e.g. Nissen & 
Schuster 2010) is the reason our initial cut in [Fe/H] for our sample 
was at [Fe/H] = -0.70 dex separating the accreted halo from the 
canonical thin disk in the [a/Fe]-metallicity diagram. Upon sepa¬ 
rating our metal-poor sample into several subpopulations in section 
2.2 we found a non negligible amount of stars that are metal-poor, 
a-poor yet the kinematics suggest they may be disk-like, which we 
called the undetermined subgroup (cyan diamonds). 

In Figure 6, it is clear that the [C-fN/Fe] and [Al/Fe] abun¬ 
dances of the undetermined subgroup at the high metallicity end 
( -0.83 < [Fe/H] < -0.70 dex) are, on average, consistent with the 
thin disk population at those metallicities while at lower metallici¬ 


ties ([Fe/H] < - 0.90 dex) the undetermined population decreases 
in both [C-bN/Fe] and [Al/Fe] and behaves like the accreted halo 
subgroup. For this reason we changed the colour of the diamonds 
with metallicity [Fe/H] > -0.83 dex to magenta in order to de¬ 
note that these stars chemically resemble the canonical thin disk 
population. In Figure 8(a), we plot [Al/Fe] as a function of metal¬ 
licity now with all a-rich stars as contours with the canonical thin 
disk (open magenta hexagons) and our candidate metal-poor thin 
disk (magenta diamonds). We leave the diamonds with metallicity 
[Fe/H] < -0.83 dex as cyan because these stars are still undeter¬ 
mined although it is likely they are either accreted halo stars or 
miscategorized thick disk stars. 

In the top panel Figure 8(b) we plot the GRV/cos(fe) as a func¬ 
tion of I and 6 as a function of I in the bottom panel for emphasizing 
the candidate metal-poor thin disk (magenta diamonds) and unde¬ 
termined (cyan diamonds) subgroups relative to the canonical thin 
disk component (open magenta hexagons). We plot all a-rich stars 
(i.e. canonical thick disk and halo) as background orange contours, 
following our notation color for the thick disk, for comparison. 
We note the following: (1) the canonical thin disk at intermediate 
metallicities (open magenta hexagons) population is predominately 
found between 90 < Z < 220 deg and (2) the candidate metal- 
poor thin disk (magenta diamonds) are spatially and kinematically 
consistent with the canonical thin disk population. Specifically, the 
GRV dispersion of the candidate metal-poor thin disk at a constant 
I (i.e. magenta diamonds in Figure 8(b)) is very small and is thus 
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consistent with thin disk. This suggests that the canonical thin disk 
may extend to lower metallicities than found in recent studies (e.g. 
Bensby, Feltzing & Oey 2014). The lower metallicity extent com¬ 
pared to other studies (such as Bensby, Feltzing & Oey 2014) could 
be a result of the fact that APOGEE has a large sampling of stars in 
the outer Galaxy, where the most metal-poor thin disk stars should 
exist. As noted in section 2.1, while the precision of the metallic¬ 
ities are high, the accuracy or metallicity scale may not be well 
calibrated. Floltzman et al. (2015) used a sample of globular clus¬ 
ters and found that in this metallicity regime ([Ee/H] = -0.6 to - 
1.0 dex), [Ee/H] can be overestimated by as much as ~ 0.10 dex. 
This causes the metallicity extent of the thin disk to extend down 
to at least < —0.83 dex if we use the uncalibrated [Ee/H] or -0.95 
dex if one calibrates the [Ee/H] scale using the globular clusters. 
We remind the reader we choose to use the uncalibrated [Ee/H] 
because the [X/Fe] trends match what is expected from literature 
while calibrating the [Fe/H] values has the effect of inflating the 
[X/Fe] (Holtzman et al. 2015). Given this uncertainty in the cali¬ 
bration, we plan to explore the metal-poor extent of the thin disk 
with other surveys, such as the Gaia-ESO survey (Gilmore et al. 
2012), in order to investigate the absolute value of the metal-poor 
tail of the thin disk metallicity distribution. 

Another interesting result to point out from Eigure 8(b) is that 
our candidate metal-poor thin disk (magenta diamonds) and most 
of the intermediate metallicity stars of the canonical thin disk, are 
found at locations outside the solar circle (90 < I < 270 deg) while 
the more metal-rich canonical thin disk is seen at Galactic longi¬ 
tudes consistent with locations inside the solar circle indicative of 
a negative radial-metallicity gradient. Negative metallicity gradient 
are widely observed in the canonical thin disk both in small sam¬ 
ples (e.g. Genovali et al. 2014) and large surveys (e.g. Cheng et al. 
2012; Mikolaitis et al. 2014). Theoretically, the presence of a neg¬ 
ative metallicity gradient in the canonical thin disk is evidence for 
an “inside out" formation scenario of the disk. Under this scenario, 
the oldest disk population is formed from the low angular momen¬ 
tum gas, which falls into the center of the halo and begins rapidly 
forming stars earlier compared to the outer disk. The rapid star for¬ 
mation leads to larger metallicities being produced in the central 
disk compared to the outer disk causing a negative metallicity gra¬ 
dient (e.g. Larson 1976; Minchev, Chiappini & Martig 2014). Our 
finding is thus suggestive of a negative radial metallicity gradient, 
which is consistent with the inside-out formation scenario. 


5 REDEFINING THE SELECTION OF GALACTIC 

COMPONENTS: A CHEMICAL TAGGING APPROACH 

The current method of decomposing the canonical Galactic com¬ 
ponent in the current scheme of chemokinematics (the use of both 
chemistry, e.g. [a/Ee] and metallicity) and kinematics (e.g. radial 
velocities, or U,V,W velocity vectors) has the disadvantage that 
it can be quite complex. The complexities lie in controlling the 
assignment of stars in the populations. Eor example, we saw in 
section 2.2 where we had to use both chemistry and kinematics 
just to separate out the Galactic components and even had an ex¬ 
tra undetermined population that could not be categorized until af¬ 
ter full chemical analysis. This caused each plot to contain many 
more symbols (and thus populations) than may be necessary. Fur¬ 
thermore, the assignment of populations on just kinematics (e.g. 
Bensby, Feltzing & Oey 2014) or both chemistry and kinematics 
causes significant biases making it difficult to discern the chemi¬ 
cal evolution of the Galaxy. This leads us to question: Is there a 


more simple approach? More specifically, can we just use chem¬ 
istry to directly separate out the canonical Galactic components to 
then study the spatial and dynamical distributions? 

In this section we attempt to answer this question by putting 
forward a chemical-only labelling approach to separate the Galac¬ 
tic components based on the above analysis. In section 3 and 4, we 
discussed that among the various elements studied here the canon¬ 
ical halo, accreted halo and thin disk were able to be separated pri¬ 
marily in the Q-elements (e.g. [a/Fe] , Mg, or O), Al, C-l-N. In this 
way, we can use a combination of these elements relative to Fe (or 
even other Fe-peak elements such as Mn) to chemically separate 
the canonical thin disk, accreted halo and thick disk. We note from 
above that the canonical thick disk and halo are effectively chemi¬ 
cally indistinguishable meaning they can, in principle, be treated as 
one population although they likely have an inhomogeneous kine¬ 
matic and spatial distribution. 

The [a/Fe] ratio is particularly important because it tracks the 
ratio of the elements produced in SNII and SNIa providing use¬ 
ful constraints on the formation timescale of the component. We 
note here that the abundance ratio [Mg, O/Mn] would actually be 
a better tracer of SNII/SNIa as Mg and O are the first elements to 
be affected as a result of SNII (e.g. Nomoto, Kobayashi & Tom- 
inaga 2013) and Mn is produced at higher fractions compared to 
Fe in SNIa (e.g. Gratton 1989). However, we choose to use [a/Fe] 
in part because it is more accessible in other non-APOGEE data 
compared to Mn as well as the fact that Mn is difficult to mea¬ 
sure in optical spectra as several lines suffer from hyperfine struc¬ 
ture splitting (e.g. Nissen et al. 2000, Jofre et al., submitted). Al 
is important as it is produced via SNII and is sensitive to the ini¬ 
tial C and N abundance, which in turn is produced by He burning 
or AGB evolution. Eor these reasons, in Figure 9 we show [a/Fe] 
as a function of [Al/Fe] in both panels. The data colored by their 
[C-l-N/Fe] in the left panel (red dots have [C-l-N/Fe] > -1-0.05 dex 
while blue dots have [C-l-N/Fe] < -1-0.05 dex) and [Mg/Mn] in the 
right panel (green dots have [Mg/Mn] > -1-0.12 dex while black dots 
have [Mg/Mn] < -1-0.12 dex). The data separates in three clusters in 
Figure 9. One being the thin disk stars with low [a/Fe], moderate 
[Al/Fe] and low [Mg/Mn] shown by a magenta selection box. The 
second is the thick disk-l-halo stars that have higher [a/Fe], higher 
[Al/Fe] and high [Mg/Mn] with an orange selection box. Lastly is 
the accreted halo stars that have low [a/Fe] with a larger spread 
than the thin disk, low [Al/Fe] and moderate to low [Mg/Mn] with 
a red selection box. 

As a consistency check of the populations defined from the 
selection boxes in Figure 9, we looked at where they fall in the 
[a/Fe]-[Fe/H] plot (Figure 10(a)). In Figure 10(a) we see that we 
recover the same structure seen in Figure 3(a). The canonical halo 
and thick disk are not chemically distinct and represent the a-rich 
stars. The thin disk is a-poor and primarily at higher metallicities 
but extend to very low metallicities and the a-poor accreted halo 
is primarily found at low metallicities although extend to as high 
as -0.60 dex consistent with Nissen & Schuster (2010). The large 
gaps in the [a/Fe] in Figure 10(a) are artificial as they are a result 
of our selection boxes. We plan to explore a more sophisticated 
statistical approache to improve or abandon the selection boxes in 
favor of a probabilistic and clustering approach. Here we show that 
in a simple case of using [a/Fe], [Al/Fe], [C-l-N/Fe] and [Mg/Mn] 
we can indeed chemically tag the bulk of the Galactic components 
in a robust way. 

This is confirmed when looking at the l-GRV space (Figure 
10(b)). We use Z-GRV space to determine if this purely chemical 
definition of the canonical Galactic components recover the kine- 
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Figure 8. (a) [Al/Fe] as a function of metallicity of the candidate metal-poor thin disk (magenta diamonds). All Q;-rich stars are shown as orange contours. 
The canonical thin disk sample is shown as open magenta hexagons for comparison. Dotted lines are placed at the metallicity cut [Fe/H] = -0.83 dex and a 
guiding line is placed at [Al/Fe] = 0 dex. (b) Upper panel: The GRV as a function of I for the candidate metal-poor thin disk. The full APOGEE sample that 
pass our quality cuts are shown as black contours. For consistency the blue dotted lines represent the initial kinematic selection. Lower Panel: The Galactic 
coordinates (/, b) of the vaiious subgroups following the same symbols as the top panel. We found that the candidate metal-poor thin disk with metallicities 
down to ~ -0.83 dex have kinematic-spatial coherence that is consistent with the canonical thin disk at these metallicities. 
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Figure 9. Left Panel: The [a/Fe] abundance ratio as a function of [Al/Fe] for all intermediate metallicity stars with —1.20 < [Fe/Fl] < -0.55 dex coloured by 
[C+N/Fe] where red is [C+N/Fe] > +0.05 dex and blue is [C+N/Fe] < +0.05 dex. Right Panel: The same as the right panel but now colour coded by [Mg/Mn] 
where green is [Mg/Mn] > +0.12 dex and black is [Mg/Mn] < +0.12 dex. In both panels, the error bar represents the typical uncertainty in the abundance 
ratios. There are three clusters within the data. The canonical thick disk + canonical halo, canonical thin disk and accreted halo components which should 
fall within the orange, magenta and red selection boxes, respectively. Note with just these four chemical dimension each component is marked by a different 
colour (e.g. accreted halo it is blue indicating a-poor, Al-poor, and [Mg/Mn] is mixed while the canonical thin disk is Q-poor, Al-rich and [Mg/Mn] is low) 
and thus will very strongly separate in these chemical spaces. 


malic structure we expect them to have. Like in Figure 3(b), the 
dotted blue line represents the sinusoidal pattern of stars following 
disk-like motion. The canonical thin disk (open magenta hexagons) 
not only have spatial-kinematic coherence indicating co-rotation 
with the disk, they are also located primarily outside of the solar cir¬ 
cle. The accreted halo (red triangles) have hot kinematics although 
some are located inside a region of Z-GRV space where they may 
be co-rotating with the disk (i.e. inside the blue dotted line in Fig¬ 
ure 10(b)). Full 3D kinematics would be very helpful to validate 
these stars as accreted halo. In addition, the accreted halo may ex¬ 
tend to metallicities as high as [Fe/H] < — 0.6 dex. The canonical 
thick disk + halo (orange stars) are dominated at these metallicities 
by the thick disk and thus by in large show a co-rotation with the 
disk. The GRV/cos(fe) dispersion is much larger than the thin disk 
as expected. However there are orange stars that reside quite far 
outside the disk-like patter (blue dotted line) which are the canon¬ 
ical halo stars. This test verifies that by using the [Q/Fe]-[Al/Fe] 
diagram in connect with [C+N/Fe] and [Mg/Mn] we can distin¬ 
guish Galactic populations in a way that can be studied free of 
biases in the dynamical spaces. This has the advantage to tradi¬ 
tional methods that full 3D kinematics are not necessary (outside 
of validation) and one can break the degeneracies in the [a/Fe]- 
metallicity diagram (i.e. the initial undetermined group is now cat¬ 
egorized, see Figure 10(a)). We also have less contamination in that 
we no longer have thin disk stars that have extreme velocities (i.e. 
the open magenta hexagons outside the dotted blue line in Figure 


3(b)). The population assignment also works just using the [a/Fe]- 
[Al/Fe] diagram alone rather than adding the addition [C+N/Fe] 
and [Mg/Mn] abundance ratios. The advantage to adding those ra¬ 
tios is to reduce contamination from the fringes of the thick disk 
and thin disk components. This is just a first step in showing that 
chemical-only approaches may be the way forward in decomposing 
the Galaxy especially with large chemical surveys either underway 
or planned for the near future. More work will be needed to move 
beyond the selection boxes we use here to more robustly quantify 
the discreetness between the components in these chemical spaces, 
which we plan to explore in a forthcoming paper. 


6 SUMMARY 

We have studied the chemistry of a relatively large sample of in¬ 
termediate metallicity ([Fe/H] < -0.55 dex) giant stars using the 
APOGEE survey in order to explore the disk-halo transition though 
a chemical tagging approach. We selected a sample of ~3200 in¬ 
termediate metallicity stars using a set of quality control cuts de¬ 
scribed in section 2.1. We used the [a/Ee]-[Fe/H] diagram to split 
and Z-GRV diagram (e.g. Majewski et al. 2012; Sheffield et al. 
2012) to split our sample into five subgroups: (1) canonical thin 
disk (open magenta hexagons), (2) canonical thick disk (orange 
stars), (3) canonical halo (black squares), (4) accreted halo (red tri¬ 
angles), (5) undetermined (cyan diamonds) as described in section 
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Figure 10. (a) The [o/Fe] as a function of metallicity for the intermediate metallicity sample using only chemistry to separate the components. The canonical 
halo + thick disk stars (orange selection box of Figure 9) are shown as orange stars. The canonical thin disk (black dots in the magenta selection box of Figure 
9) are shown as open magenta hexagons. The accreted halo (blue dots in the red selection box of Figure 9) are shown as red triangles, (b) Upper panel: The 
GRV as a function of I for the intermediate metallicity sample using only chemistry to separate the components with the same symbols as (a). For consistency 
the blue dotted lines represent the initial kinematic selection. Lower Panel: The Galactic coordinates (Z, b) of the various subgroups. 


2.2. In section 3 we discussed the chemical abundance patterns of 
these four groups in a-elements (Mg, Ti, Si, Ca, O, S), Fe-peak el¬ 
ements (Mn, Ni), and odd, even, and light elements (C-l-N, Na, Al, 
K, V). 

Our main results of this Galactic chemical evolution analysis 
can be summarized in the following points: 

• The canonical halo and canonical thick disk are not chemically 


distinct in any of the elements studied. This is most evident in the 
abundance ratio distribution in the various elements shown in Fig¬ 
ure 7 which indicates that the a-rich halo component (black line) 
and the a-rich disk component (orange line) overlap both in terms 
of mean and dispersions for both populations. This may suggest 
that the Galactic halo and thick disk are chemically formed from 
similar gas but one is pressure supported and the other is angular 
momentum support. 
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• We uncovered the largest sample of ‘accreted’ halo stars. We 
showed, for the first time, that accreted halo stars are under abun¬ 
dant in [S/Fe], which is unsurprising given its status as an a- 
element. We have also showed that the accreted halo is under abun¬ 
dant in [Al/Fe] which is likely caused by the nuclear pathways in 
which A1 is produced. The low abundances of [C-l-N/Fe] of accreted 
halo stars respect to canonical stars imply that this gas had a slower 
chemical environment (Nissen et al. 2014). The Ni of the a-poor 
sequence is, on average, under abundant relative to the a-rich se¬ 
quence consistent with previous studies. These results indicate that 
we have confirmed the existence of an a-poor halo sequence of 
stars that is chemically different than the stars formed the rest of 
the halo with chemistry suggesting that these stars may bom ex¬ 
ternally (e.g. Nissen & Schuster 2010, 2011; Schuster et al. 2012; 
Ramirez, Melendez & Chaname 2012; Hawkins et al. 2014). 

• Upon closer investigation of the undetermined population, we 
have found evidence that the thin disk extends down to very low 
metallicities ([Fe/H] < —0.83 dex). In Figure 8(b), we showed that 
the candidate metal-poor thin disk, which was the metal-rich end 
of the undetermined population, were consistent with the canonical 
thin disk at comparable metallicities both kinematically and spa¬ 
tially. This was the first indication that these stars may be among 
the lowest metallicity thin disk stars. However, this result will need 
to be confirmed within another survey with a more accurate metal¬ 
licity scale. Spatially, the candidate metal-poor thin disk as well as 
the intermediate metallicity thin disk is located in the outer Galaxy 
compared to the more metal-rich thin disk component located in the 
inner Galaxy consistent with a negative metallicity gradient which 
is observed in the literature (e.g. Cheng et al. 2012). 

• Finally, in section 5 we put forward a powerful chemically 
labeling approach to separate the chemically distinct canonical 
thin disk, thick disk-l-halo, and accreted halo components using 
the [a/Fe], [Al/Fe], [C-l-N/Fe], and [Mg/Mn] abundance ratios. 
For consistency in Figure 10(a) we plot our separated components 
on the standard [a/Fe]-metallicity plane. Finally, we verified this 
method by confirming that the kinematics of the resulting popu¬ 
lations matched the expected distribution in Z-GRV space (Figure 
10(b)). 

Ultimately the precise proper motions and distances from Gaia 
will allow us to study the full spatial position and dynamics of 
the stars from each of the components to compliment the chemi¬ 
cal evolution analysis done here and help solidify the existence of 
the accreted halo and candidate metal-poor thin disk. Future sur¬ 
veys such as 4MOST (de Jong et al. 2012) and ongoing surveys 
such as the Gaia-ESO survey (Gilmore et al. 2012) and GALAH 
(De Silva et al. 2015) will be able to further study the chemical pe¬ 
culiarities of the accreted halo in a statistical way. With full space 
motions, positions and chemistry of a large sample of stars in hand, 
we may not only be able to study the accreted a-poor sequence 
but also confirm the existence of the very metal-poor thin disk and 
chemical continuity between the halo and thick disk. 
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